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Abstract 

We study the microcanonical description of string gases in the presence of D-branes. 
We obtain exact expressions for the single string density of states and draw the regions 
in phase space where asymptotic approximations are valid. We are able to describe the 
whole range of energies including the SYM phase of the D-branes and we remark the 
importance of the infrared cut-off used in the high energy approximations. With the 
complete expression we can obtain the density of states of the multiple string gas and 
study its thermal properties, showing that the Hagedorn temperature is maximum for 
every system and there is never a phase transition whenever there is thermal contact 
among the strings attached to different D-branes. 



1 E-mail address: laucelli, jesus@stringl.ciencias.uniovi.es 



1 Introduction 



Finite temperature frameworks are in general very interesting laboratories in which to 
study the fundamental degrees of freedom of a theory. In the case of string theory the 
interest increases because of the exponential growth of the number of states with the 
mass. This behavior, known as the Hagedorn spectrum is relevant because of the 
fact that it generates a critical thermodynamics that in principle could correspond to a 
phase transition [0, [31, |l H. The analysis of this possibility was firstly attempted from 
the canonical ensemble [0, [3], [4], |J, but it seems that a more fundamental description of 
the degrees of the system is needed. 

The description that, in principle, would clarify the open questions on the Hagedorn 
thermodynamics is that of the microcanonical ensemble J?], [| §, || . The density of states 
of the string gas was studied in the early works [[I], || on the subject by taking an am- 
biguous high energy approximation. In fact the thermodynamics for the closed strings 
that was derived there do not have any information about the low energy behavior of the 
system. This prevents any trustable study on the way the system reaches the Hagedorn 
temperature. 

On the other hand in [0, |6| there was a disregarding of the volume dependence on the 
density of states that led to the surprising conclusion that there is a difference between 
a system with an already open universe and one in which we take to infinity the radii of 
the initial torus 0. 

In |8], another point of view was claimed. The density of states of the single closed 
string valid for all energies was obtained. The multiple string one was formally defined 
from the convolution theorem || , and it was shown that for an already open universe the 
systems reaches the Hagedorn temperature for a finite energy per string. At this point 
a phase transition occurs in such a way that the temperature is kept constant and the 
specific heat diverges. By this analysis it was possible to know how the systems evolve 
from low energy to the high energy limit. 

The compactified case was studied in |J . In this background the system has a unique 
phase with finite and positive specific case. The Hagedorn temperature works as the 
limiting one. An interesting point was that in the decompactification limit the system 
coincided with that of |§. 

A common feature appearing in both cases is that at high energies equipartition breaks 
at a given energy per string. From this point the system is composed of a dominant highly 
energetic long string in thermal contact with a sea of a large number of low energetic 
strings that works as a thermal bath. This relevant part of the system was not taken 
into account previously in |7], || because it emerges from the IR part of the single string 
density of states. The existence of the string sea is relevant in order to prevent the negative 
specific heat phases. 

From an completely different perspective the same picture has been recently obtained 
fl6f . The authors of |l(J, considered the possibility of a high energetic radiating string. 
They discovered that this string would behave as a black body at a temperature T H and 
if one included a thermal bath composed by low energetic strings the system would reach 
equilibrium only if both temperatures, those of the black body and those of the bath, 
coincide at the Hagedorn value. This strongly supports the idea presented in §, §], that 
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a long energetic string would be in thermal equilibrium only in the presence of a thermal 
bath at T H . 

To put it in short: working with approximate expressions without a neatly defined 
range of validity is dangerous. It is difficult to discern in which situations low energy 
effects can qualitatively change the conclusion of the physical analysis. This is specially 
true in systems composed of many objects, let alone those where a thermodynamical limit 
is to be taken and even more near a phase transition, because local and small fluctuations 
compared to the total energy can always make individual objects visit the regions of phase 
space that have been thrown away. In our work we will render an exact calculation of 
the single string density of states and a comparison with the different approximation that 
come from it will be done. 

These considerations hold for perturbative string physics. It has been sometimes 
argued that non-perturbative objects would have some role to play in relation to the 
Hagedorn problem. There have been some works in which this topic has been dealt 
with, either trying to include D-branes |15|] in the canonical ensemble [II], |I2| , or 



D-instantons [T3| effects in the microcanonical ensemble. A laconic 'nothing happens' has 
been the usual conclusion. In the former case the problem of approaching the Hagedorn 
problem from the canonical ensemble is exactly the same that for closed strings. That is 
the canonical description breaks and nothing else could be said since this point. In the 
latter the main argument is that the energy scales at which the non-perturbative effects 
become important are completely dominated by the perturbative Hagedorn physics, so it 
seems natural to conclude that a subsequent effect cannot regularize an already existing 
at low energies one. 

D-brane dynamics is governed by the Dirichlet open strings attached to it [pl|, so in 
order to study D-brane thermodynamics we need to study the finite temperature version 
of these open string theories |K| [U], [I2| . On the other hand if we want to study the 
complete thermal behavior of type II string theories we need to account for the closed 
string sector and for the corresponding D-branes allowed for each theory. Finally if we 
want a fundamental understanding of the Hagedorn 'phase' of the theory we need to study 
it in the microcanonical ensemble. This idea has been recently developed in 



In this work the microcanonical density of states has been computed for a variety of 
brane configurations by using the methods of |6|, [7|], and high energy thermodynamics has 



been described. A related work was also made in [14 . 

In the present paper we approach the microcanonical description of thermodynamics 
with D-branes as background applying the methods of || [|. The plan of the work is 
the following. We begin the analysis by computing the density of states from the mass 
spectrum. The next step consists of a computation of the single string density of states by 
the method of the inverse Laplace transform. We separate the analysis for high energies, 
analogous to the analysis of fTjjfl, from those valid for the whole energy range. Using the 



previous results we go to the thermodynamical description of the system . 



2 The Density of States from the Mass Spectrum 

There are two ways to compute the density of states of the single string system. The 
first one consists of the application of the inverse Laplace transform. The method begins 
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with the computation of the Helmholtz free energy and then using it as the single string 
partition function in such way that we obtain the density of state from the Bromwich 
integral 

nW(£) = / dpzw)*?* (1) 

im Jc 

where Z\(f3) = —Vj3F(/3). The n-string term is characterized by the partition function 

Z n ((3) = -Mf3) n . (2) 

The other way to calculate the microcanonical density of states of a system is to 
count the degeneracy in the phase space geometrically. It is specially useful in the case of 
compact target backgrounds. In this section we shall implement this idea in our case of a 
gas of open superstrings moving in the presence of D-branes. We start with the following 
Hamiltonian, that corresponds to the mass formula of the theory which. We take it as a 
dispersion relation 



a'E 2 = p 2 + u 2 + 2N (3) 

where p is the momentum, that is discrete in this case. The dimension of this vector is 
the number of Neumann directions, Q is the winding number. It is always perpendicular 
to the momentum vector because open strings cannot wind around a direction and have 
a momentum at the same time. The dimension of this vector is, thus, the number of 
Dirichlet directions. N is the total oscillation number. It is the sum 



N = J2nN n . (4) 

This means that there are several states with the same total oscillator number. They 
are counted by the following coefficients: 

~~ 2rri f q k+l \T] 3 (q) J {^(q) ) [ ) 

where the integrand is the internal partition function of the theory 

Let us begin to calculate the degeneracy associated to a value of the energy fixing 
firstly iV to a constant. In that case: 



On the other hand 



f + u 2 = a'E 2 — 2N = C (6) 



/ + ^=E^+Ef (7) 



2 The dependence of the internal partition function of the theory on the number of mixed, ND and 
DN, direction comes from the different normal ordering constant of these sectors. We are very grateful 
to J.L.F. Barbon for focusing our attention on this point. 
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where NN marks the directions with Neumann boundary conditions and DD the Dirichlet 
ones, there are also mixed directions, DN + ND = u, that do not contribute to the kine- 
matical energy but through their oscillator masses only. This contribution is accounted 
by the coefficients a at . To make things more homogeneous, we can use the variables 



if % is Neumann 



so that 



\fa' 

Ti = —j— if i is Dirichlet 



this is the equation of an ellipsoid in the discrete phase space. We are interested in 
counting the number of states for which that relation approximately holds, that is, letting 
the energy be in [E — AE, E + AE]. If the energies are very low, it is important to remind 
that AE has a lower bound which is the maximum gap between two energy levels. The 
number of states is proportional to the volume of the ellipsoidal shell defined by equation 
(P). The total volume is 

v(C) = n 3 (j[r?\c £) /*. (10) 

The volume of a thin shell is 

= VL b D \J[ rA (a'E 2 - 2k) ^ a'EAE (11) 
where Qjj is the angular part. Its dependence on the dimension is 

n n = _ — = — (12) 

D DT(D/2) V ; 

In our case D = DD + NN < 9 — u, is the number of testable directions. 

To get the total degeneracy it is necessary to multiply by the one associated to iV (ajv) 
and by the size of the super-multiplet: 2 4 . After the sum over all possible values of iV is 
performed, the result is 



fl d (V; E) = 2 4 ^L (2^)^"^^ a ^'E (a'E 2 - 2N) 2 H (a'E 2 - 2N) 

r(D/2) Vdd n=0 ' 

(13) 

This is not valid for D = because the derivative in (|TT|) is not correct. In that case 
there is no continuous contribution from either the momenta or the windings and do the 
density of states is just a sum of delta functions centered on the mass levels. We shall 
talk further in section 5. 



DC 
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For a value of E, the sum over N is finite because the oscillator number is limited by 
the available energy. This is expressed by the Heaviside function. The calculation will 
be complete when we take into account the quantum statistics. Before doing that, let us 
add some terms that are important if the proportion between some radii is very different 
from one. In that case, there is a range of energies where the first non-zero modes of some 
directions lie far beyond the ellipsoid in the phase space and the geometric count is no 
longer valid because the single zero mode is not counted as one but as a small fractional 
number of states. To simplify, we shall consider all the Neumann radii equal and the 
same with the Dirichlet ones. The only parameters are then Rdd and Rnn- This can be 
generalized to consider all the possibilities but there are too many and all the interesting 
physical consequences can be extracted studying this simpler case. The density of states 
that this yields is 



4 2n D/2 /r ,, DD _ NN V NN 



oo 



fLi (V, E) = 2 4 —„ (2ir) DD - NN ^ V a N a'E 

cl rp/2) Vdd^o 



-1/C 

'c 



x (a'E 2 - 2N) 2 H (a'E 2 - 2N) H(E-V c 

9 C/2 (n^-\C oo c_i 

+2A wrj^-^- £ aNa ' E ( a ' E2 - 2N ) T ~ n ( a ' E2 - m ) n fa - E 



r(c/2) v c N=Q 

where 



(14) 



If R NN ^>1/R DD C stands for NN and V c = V NN 

If R NN « l/R DD C stands for DD and V c = l/V DD (15) 

Now it is easy to add quantum statistics following this recipe 

-(3F((3,V) B , F = ±j:iog(lTe^^) = 

{n} 

= T ^g(±ir^ ({ n } ) (16) 



{n} r=l r 



This leads to the final result 



9 D/2 T/ oo oo I 

n(V,E) = 2^^L— ^^nnVnn^ E J- aN a'Ex 



— — 1 

x (a'E 2 - 2Nr 2 ) 2 H (a'E 2 - 2Nr 2 ) + 

n C/2 lc\ \C oo c i 

+2 \(C/2) Vc f o a Na'E(a'E 2 -2 N y- H (a'E 2 - 2N) H (V c 1/C - E) . (17) 

This result is completely analogous to that of 0. We will see in next sections how we 
can recover the same density of states for the single string system from the study of the 
free energy of the theory in the canonical ensemble. 
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3 High Energy Density of States 



In this section we are going to compute the Helmholtz free energy of the Dirichlet open 
string theory for a variety of D-brane backgrounds from the canonical ensemble, and we 
will finally study the different limits of the density of states derived from it. We will do 
it for the asymptotic high energy range. The approximation is based on the idea that 
the existence of the Hagedorn critical temperature is an UV phenomenon so it should be 
enough to look at the very excited string states to obtain the relevant information, this is 
the perspective assumed in H, fF|, and more recently, for the D-brane backgrounds, in 



The other option is to try to compute the complete density of states without projecting 
out any string states, and will be detailed in the section 5. 

We will begin by noticing that the degrees of freedom lost in this high-energy limit, 
are, in fact, projected out in the canonical description by a IR cut-off that cannot be 
forgotten. If we introduce it explicitly, we would easily see how the massless open strings 
are neglected. This projection prevents the possibility of accounting for the low energetic 
sea. 

Let us firstly fix our background. Suppose we have a configuration with a Dp- 
brane and a Dq-brane. We can define four different types of directions, the Dirichlet- 
Dirichlet(DD), Neumann- Neumann (NN), ND and DN. The last two do not contribute to 
the mass of the strings, apart from the oscillation terms (That means that we can neither 
have momentum nor winding in these directions). We will include in this classification of 
the directions the spatial ones only. They are related by 

DD + NN + DN + ND = 9. (18) 

The Euclidean time is supposed to be compactified on a circle of length (5. Another 
interesting parameter we will deal with is the sum of the mixed directions, v = ND + DN. 
The interest of it comes from the fact that the when v = 0,4,8 the string gas is evolving 
in a completely supersymmetric background. We will focus our study on this type of 
brane configuration, however some comments on more general ones will be made. 

In brief, the string's quantum numbers are: the momentum mode in the NN directions; 
the winding mode in the DD ones and the Matsubara frequency in the Euclidean time, 
the DN and ND direction do not contribute to the mass of the strings. 

We can now compute the Helmholtz free energy as the vacuum energy of the theory 
described above. It formally reads 

F f°° dt 2ira'tM 2 



f3F(f3) = -Tr H (-lY e -2natM^ (ig) 

J It 



where M is the mass of the string and is given by 

NN 2 DD A-K 2 m 2 
M 2 = T.% + £«M + ^ + Osc (20) 
i=i K i i=i P 

where the last contribution comes from the oscillator modes of the open string. We can 
compute the trace over the Hilbert space obtaining 

- pm = -J f n 0a k n * a (o, ^r\) m (21) 



where the internal partition function f(t) is given by 



The high energy behavior of the open string gas may be studied in the canonical 
ensemble and that is the easiest way of visualizing the Hagedorn critical behavior. Let us 
begin our study by looking at this region. When the energy is large the string modes that 
dominate the dynamics will be the massive ones. These are enclosed in the UV region of 
the previous expression, that is accounted by the large t limit. We can introduce an IR 
cut-off A 2 in ([21]), and in this way we will forget about low energetic strings, we have 



A 2 df nn / oita'\ DD / \ ( Rir 2 ita'\ 

-w)=-jf In^^^jn^^^ito^^^^j/w (23) 

In order to study the UV limit (t —>■ 0) it is convenient to take advantage of modular 
properties of the theta functions in the free energy going to the closed-string channel by 
making t = 2n 2 s~ 1 , obtaining 

V NN R r°° NN / ;d2 \ DD / • \ 

w> ~ -hnl L * n * (o, ^) n h (o, ^) 





POO 


'Air. 








x6 4 





We can now go to the large s limit by simply taking into account that 



(24) 



s^oo 



we finally can express the UV limit of the open string free energy, with the correct mul- 
tiplicative prefactors, as 



\/NN a roo 

mP) = -^W DD - NN W)^£-(2nr i E J 2n2 dse~^ M ^ (26) 



1 ,1)1)- XX,,, i!) ' ( -,— j ■ 

{M,n} J ~ 

where we have defined 



n 2 B 2 1 

A(M,n,/?)=M 2 + — i-,-—, (27) 

in that way the winding and momentum modes are included into the mass M. In what 
follows we will collect the numerical factors into a multiplicative constant C. When 
A(M,n, (3) = we arrive at a critical temperature analogous to the Hagedorn one. In 
fact we can see that at those temperatures a string state becomes massless. When the 
temperatures are passed through this states get a negative squared mass, that is, it 
becomes tachyonic. We can solve the equation for the critical temperatures obtaining 

„ 9 Sua' 16ttVM 2 , n 

(3 2 = — with < n G Z. (28) 

n 2 n Z 



8 



With n = 1 and M = we have the lowest critical temperature that corresponds to 
the Hagedorn one. The following critical points are obtained for bigger temperatures, 
in fact we can in principle have negative unphysical solutions for 3 C . This possibility is 
avoided by assuming that the compactification radii are bigger than the selfdual ones. 

At this moment we can write the Helmholtz free energy as a sum over functions which 
are analytical in the whole complex-/? plane except in the neighborhood of 3 C . The sum 
reads 



yNN a f00 

-«=-C^l/>'^'. (») 

{/3 c ,n} a2~ 

Finally we can integrate the previous expression in order to obtain a closed expression 
for the free energy in the high energy limit. Using the integral representation of the 
incomplete Gamma function we have 

~n 2 (3 2 - 8 2 c 



yNN q 

{Pc,n} 



2n\ (3 2 H 



1 ' 7TU 2 (3 2 -3 



The problem of obtaining the single string density of states from the free energy below 
is reduced to an inverse Laplace transformation. It is known that the result of this kind 
of computation is closely related with the analytical properties of the function to be 
transformed. We can look at (|50D as a sum of independent terms seeing that each one is 



characterized by its behavior in the neighborhood of its thermal critical point 3 C . We see 
that the analytical properties are enclosed in the simple pole and in the more complicate 
properties of the incomplete Gamma function. The computation of the density of states 
here would consist of a straightforward adaptation of those in || . 

Before showing the results of the computation let us comment some words about the 
cut-off A 2 . We can try to extend the validity of the approximation taking the value of A 2 
to infinity. It will give us the complete version of the Gamma function. The analytical 
properties of the Free Energy are enclosed inside the single pole appearing at each critical 
temperature. If we complete the IR part of the Free Energy by simply taking A 2 to 
infinity we will obtain an IR regular series for the density of states of the single string. 
The thermodynamics derived from it would have a positive and finite specific heat phase 
only, in which the asymptotic maximum temperature is the Hagedorn one. We think that 
this corresponds to a misconception. The IR behavior that one could obtain from the 
previous procedure is at least not complete. It is due to the fact that if we took a large 
A 2 limit we could not truncate the string theory partition function as we did in (|25|) , but 
we should have included all the string modes. 

We think that the correct way of using the previous free energies consist of keeping 
A 2 finite, and making inverse Laplace transformation of the Laurent series around the 
critical points. Doing it we will see that we obtain an asymptotic series that does not 
work for the IR regime. 

Following the previous prescriptions we can make the same computation of the high 
energy single string density of states as in ||, it gives 

yNN 1 oo (_-\\rn+ P p c E / 2tt/?T7 2 \ 

V {ft , n)F 0m=0 ^ V Ph A J 
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where A((3 c ,n,m,p) is a constant independent of E but dependent, through /3 C , of the 
R parameters, and TC(x) is the Heaviside function. The sum in m comes from a Taylor 
expansion of the incomplete Gamma function. 

The behavior is in general governed by a leading term that corresponds to the Hagedorn 
singularity. We will come back to the details later but let us sketch the general dynamics. 
We find a positive, but infinite, specific heat phase. The leading term does not give 
any negative powers of the energy taking us to theory with a constant temperature and 
infinite specific heat behavior. The subleading terms give corrections for small energies. 
The string gas takes a finite and positive specific heat phase. 

We think that these corrections are meaningless because the density of states is only 
valid for a very high energy, and in order to know the low energetic phase we need an 
IR-complete density of states. We will come back on these features later when studying 
the thermodynamical behavior of the gas. 

It seems interesting to study the high energy density of states when we open some of 
the DD directions. This case physically means that there would be a reduced number 
of testable directions for the single string and has shown its interest in previous studies 
|7|, [18], [D]]. The result is analogous to ([31]). It is implemented by taking some of the 
DD, (let it be a number < a < DD of them), infinite. We will also suppose that the 
momenta on the NN directions are continuous so 

nOO -qDD tdNN r>NN _^ „ /oo\ 

ti x ,...,lt a ,K 1 , ...,lt NN — ► OO. [6Z) 



Including it into (|2~4|), and then also in(p8"D, we arrive at the following UV canonical Free 
Energy 



yNN n 

{Pen} 



n 2 ((3 2 - (3D 



(9-D-v) 
2 



•1 



r i + 



D + u-9 Tin 2 {(3 2 -(3 2 c y 
2 ' A 2 (3 2 



H 



27T P 2 H 

(33) 

where D is defined as a effective number of testable dimensions. It is given by 

D = (DD - a) + NN <9-u. (34) 

The meaning of this number of dimensions and the existence of 'testable' energy scales 
will be discussed later. 

In this case the ananlytical properties are more complicate. In general we can see that 
we have a branch point located at each critical temperature, concretely it will happens 
for an odd number of untestable DD directions, while for even ones we would have high 
order pole singularities. The IR continuation, by taking A 2 to infinity, is not well defined 
here. The density of states that results from inverse Laplace transform (|33D is 



yNN i oo (_-\)m+ P p c E / 2irS n 2 \ 

n^*) A(A,„, m , P ) ^!>_ B) ; +m+ H (e - . (35) 

We will see in the next section that these single string systems have a negative specific 
heat phase only. 
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The marginal case corresponds to opening all the DD direction obtaining the density 
of states of a pure NN theory 



BO) (E) = V»« •£. E E W» m, P) E l_J_ NN y, Hm+ H (e - ?5%- . (36) 



{^c,n}P=0 m=0 



This result describes a purely Neumann open string theory in iViV target dimensions. 
The energy dependence of the thermodynamics in the microcanonical ensemble exactly 



matches with a certain brane configuration described by (J35) , but without any volume 
dependence of the DD side. 

In fl3l|) and following, we can see that the density of states for the single string system 
obtained with this procedure is useless in the IR regime because of the use of the IR cut- 
off. This property prevents the use of this magnitude for the computation of the multiple 
string system by the convolution procedure |5). 



4 High Energies Thermodynamics 

In this part of our work we will go on the thermal properties of the single string system by 
using the density of states in fl35l) . We will begin with the analysis of the ranges of validity 
and we finally try to extract from it some results on the multiple string thermodynamics. 

In this section we will be interested in comparing our results with those of [18| and we 
will find where their asymptotic approach breaks down making the convolution approach 
of H, H indispensable. 

Let us begin by simply computing the microcanonical temperature. At first order it 
reads as 

> - ~ h H - H * i E - a ■ <*> 

Where we have assumed that the dominance of the term corresponding to the Hagedorn 
singularity. The properties that we can derive from it hold only at high energies. We will 
detail the analysis of the validity range in detail later. 

The temperature that results from ( p?T| ) are drawn in the Fig.(l), where we have plotted 
the temperature versus the energy of the string. We can see that the systems have negative 
specific heat phase in all cases except for D + v = 9 in which the entropy is proportional, 
at leading order, to the energy. In this case the system has an infinite specific heat phase 
only, that corresponds to a constant temperature. This is the only case in which the 
Hagedorn temperature plays the role of limiting temperature for the single string. The 
extrapolation for the multiple string system must be done carefully^. We have introduced 
the effective dimension D that correspond to the number of dimensions the system can 



In order to connect with the notation of [11 811 remember that 



7=^^-1. (38) 
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Figure 1: Temperature versus energy of the various single string system. We see that we have 
positive specific heat for the D + v = 9 case only, that corresponds to a system that tests the 
whole target space. 



probe. In fact we have the NN dimensions, that will be supposed to give a continuous 
momentum spectrum, and DD — a Dirichlet dimensions the string can wrap around. To 
eliminate possible winding excitations around a Dirichlet directions we have already taken 
the corresponding radii to infinity. However we could wonder that the system has not 
enough energy to reach a winding state even when the radii are finite. Physically this 
situation can be clearly explained for the single string system but it comes to us more 
subtle for the multiple string gas. In fact in order to prevent the winding excitations one 
has to take the energy lower than R^. If we impose it for the multiple string gas it would 
reduce the system to evolve in a very low energetic phase. 

Let us see it explicitly. If we introduce the parameters pn, that corresponds to the 
energy per string, and the number of strings per unit of NN volume n, we conclude that 

R 

pN < = puJ - (39) 

Then if we want to take a thermodynamic limit (V NN — * oo, and all the densities finite), 
we need to take the a radii to infinity too. 

In summary, if we want to make a dimension untestable for the multiple string system 
we should to take its radius to infinity, if we did not do it we might be able to talk 
about reduced effective dimension at very low energy densities p^. In order to make our 
work complete we will make the analysis for string gases with large, but finite, R a radius 
specially focusing our attention on the its limitation. 

Firstly we need to find the energy ranges in which the term corresponding to the 
Hagedorn critical point in fl35|) dominates the thermal properties. It will correspond to 
those energies at which the subleading terms are negligible, that is 

e -V> a -f)c)E „ ( O ) ^ {J3h - Pc)E > 1. (40) 
On the other hand we need to remember that the Hagedorn term in ( |35|) is null for 

where this condition comes from the Heaviside function. The combination of the previous 
expression gives the range of energies at which the system has a simply Hagedorn behavior. 
In order to specify it we need to know how is the dependence of fl4"0|) on (3 C - 
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Figure 2: The phase space of the single string system that comes from the Hagedorn dominance 
rule ( f40| ) and the winding excitation energy E w 



We know that the critical points are given by (]28|) . Because of the assumption that 
momenta should be continuous we only deal with critical points coming from winding 
modes and high-mass levels of the string. For the latter we conclude that the second 
critical point fa = \phwt' could be neglected at energies 



so we must worry about critical temperatures lower than (3% only. The first winding 
critical temperatures corresponds to the first excited winding mode and its influence on 
the thermal behavior is suppressed at energies 



Finally we will also consider the energies at which the winding states become to appear 
and consequently the density of states that neglects it becomes to fail. It occurs when 
E > Rdd/oi' = E u . Along all these consideration we have assumed the radii of all the 
Dirichlet directions to have the same value Rdd- 

We are now able to consider the phase space defined by Rdd and E and to study 
the different regions that emerge from the previous conditions. These regions will be 
separated by the equality relationship between E u , Eh and the energy of the system E. 
The resulting phase space is sketched in Fig. (2). 

Suppose the system we are interested has DD — a Dirichlet dimensions compactified 
at radius Rdd-, while the purely Neumann and a Dirichlet ones are already open. 

If Eh <C Eu <C E the winding states are excited when the Hagedorn term of the 
density of states clearly dominates the dynamics of the system. As a consequence the 
winding influence could be neglected. In this case the thermal behavior of the single 
string is accounted by the first term of ([35|) with D = NN + DD — a. This is in fact the 
situation in which the radius of the DD direction is small and could be better treated by 
tacking advantage of T-duality converting the DD directions into NN ones. In the figure 
this regime correspond the region /. 



E > 




Eh- 



(43) 



v/^ra 7 (ya'R DD - (a'R 2 DD - 1) 1/2 ) 
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In the case in which Eh <C E <C E u , the region II of the picture, we have a system 
that effectively lives in a space without the DD directions and whose thermodynamics is 
governed by the Hagedorn term. The density of states in this case is the leading term of 
(|36|) . However the existence of this region is closely related to the value of the cut-off A 2 , 
and appears at energies E ~ E c , so we think that it must be studied by the IR regular 
density of states. 

We can go now to the region /// in which E <C Eh <C E u . This region behaves 
exactly in the same way as the region //, but the difference comes from the fact that it 
appears for every value of the cut-off A 2 . On the other hand it must be described by the 
high energy density of states in an a-dimensional background. 

Finally we must look to the region IV. It corresponds to E u <C E <C Eh- In this 
case all the terms of (|35|) are relevant. We think that it would be better accounted 
with a complete density of states, overall when going to the multiple string system. The 
advantage of the asymptotic expression over the exact one, that is simplicity, is lost and 
it is better to use the complete series. 

We would like to make two final comments on the thermodynamics at high energies. 
It is interesting to note that if we pass from the region /// to the I of the phase space, 
by increasing the energy of the system, it could test large extra dimensions. Then by this 
mechanism we can change the dimensionality if the background in which the system is 
evolving. 

On the other hand it seem remarkable to us that there are some target background 
in which by this transition, we pass from a system which has a negative specific heat to 
another one with a positive and finite specific heat and which asymptotically ends at the 



Hagedorn temperatures. In terms the nomenclature of |18| we go from a non-limiting case 
to a limiting one. 

Finally we would like to mention how these densities of states could apply for the 
multiple string gas. Using the convolution description of || we know that the multiple 
string gas density of states is obtained as 

OO 1 

n(E)=52—nM(E)*.-.*nM(E) (44) 

N=0 iV • 

where the stars stand for convolution integrals. One of the possible approximations at 
high energy that we can do for the density of states of the iV-string gas is to consider it 
as behaving as a single long string. If we use this picture we arrive at 

1 1 / v NN \ N 

M na)(£0*...*nW(£?)__^^J sW{E) (45) 

where A is a factor that must be computed for each case. The simplest one are those in 
which the single string is Hagedorn dominated, and we finally obtain an ultra-high energy 
density of states as 



V NN 



It corresponds to the leading order in the computation of |7|, [18| and the corrections come 
from the series of fl^(E). The thermal behavior of the multiple string system exactly 
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matches with the single string one with the previous assumption, as for fli(E) it does 
not work in the IR region. We finally conclude that only a complete computation of the 
single string density of states can be trusted for the multiple object gas. 



5 The Complete Density of States 

In what follows we will obtain the complete density of state for the single string and use 
it to describe the string gas thermodynamics. The idea here is to try to inverse Laplace 
transform the whole free energy of the open string gas. The method that we will follow 
was presented in ||. We write the free energy of the open string theory as 

rNN 



- (3F((3) = ^2aT {NN+DD+m ^— E / -r( NN+DD+i y 2 e-% M2 f(t). (47) 

V 4?r {M} J ° * 



where the internal partition function is taken as in (|22|) . The parameter M 2 is the mass 



of the closed string boundary states propagating between the branes with the Euclidean 
time compactified on a circle of length f3. It is given by 



M 2 = - 

a! 



■NN flfl m 2 

R 2 ' 4vr 



E^ + E^ + ^ + i/2) 2 



(48) 



At this moment we need to express the internal partition function as a sum over open 
string states. More concretely we want to explicitly write 

f{t) = f j a k e-< t with a h =-±-J-% £ f(q). (49) 

fc=o Z7U J q 

Including the previous Taylor expansion and the parameter v in (fT7D, we can finally write 



PF{0) = H2aT {NN+DD+1)/2 ^mJY- EE«fc / ^r^/Vi^V^. (so) 

V {M} k=0 J ° 1 



By using the integral representation of the Bessel functions K(z) we can express the 
Helmholtz Free Energy of the open string gas as a sum over the Free Energies of each 
target massive fields in the following way 

*mJ- , (51) 
v ^ J 

where the constant in front of the integral (fi?D have been ignored. 

At this moment we would like to obtain the single string density of states by the inverse 



Laplace transformation of the Free Energy in (51) as it was done in ||. The complicate 
dependence of the argument of the Bessel function on the compactification parameters, 
the radii R, and on the inverse temperature (3, makes this computation hard. We will 
attempt it by a simple trick, that however seems useful to study some interesting and 
simpler limit before. 



15 



Let us begin with the study of the decompactification limit that corresponds to tak- 
ing all the radii to infinity. In this background the relevant degrees of freedom would 
correspond to the momenta in the NN directions. As we discussed in section 3 this is 
the only way to completely avoid the winding contribution. It is due to the fact that we 
can forget about the winding states for the single string by taking its energy smaller than 
the squared Dirichlet radius Rdd, but in order to do the same for the multiple string gas 
we must prevent the statistical fluctuations that could allow a winding nucleation. The 
density of energy is finally constrained to be in a low range of values by ([39]). In summary, 
the correct way to study the system in absence of winding states consists of taking the 
radius of the DD direction to infinity and not to assume string energies lower than the 
winding masses. 

We take the desired limit in (pi] ), by assuming the NN momenta to have a continuous 
spectrum and taking only into account the winding zero mode we arrive at 



- /3F(/3) ~ ~V NN ^ £ £ a^kn^K^ (nM n V2k) ( ttM^ ] , (52) 

where the masses now are reduced to 




iVJV+1 
2 



M„ = (53) 



'a'n 



The density of states corresponding to this free energy can be obtained by using exactly 
the same method of M, getting the following result 



n< >(E) = r( ™, EE<*. « [ E ~ ^) ■ <«) 

We can compare the high energy density of states for the single string system in the 
open universe case, that corresponds to the result expressed in (^) . The crucial difference 
comes from the IR behavior. The expression in ([54]) is regular at low energies because we 
did not exclude out any mass levels of the string, as we did for the first case, by using A 2 . 
The critical Hagedorn phenomenon here is hidden in the growth of the coefficient a&, and 
appears as the energy is big enough to open massive string thresholds. 

We can compute the density of states for the string gas living in a universe with 
compact dimensions. The way of computing it is the same we have used above. In this 
case the free energy is represented by the K\i<i(z) Bessel function. The trick that we will 
use here to use the result of || and obtain the density of states consists of writing the 
free energy in ( f£7D as 



-f3F{f3) = -{2 a 'Y {NN+DD+1)/2 — Y Y n-t- 1/2 e-£-* in+1/2)2 e-% M L w (55) 



13 £ y fC ° dt - 

™=0 {M knw } ■ 

where the mass threshold here are expressed in terms from the open string side as 



NN a 2 DD 2 

M ln W = E"S- + E H R ) + 2^7 k - (56) 

i=l n i j=l a 
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The case in which all the dimensions are compact and their radii are small enough to 
force us to consider the spectrum of momentum and winding strongly discretized, is very 
subtle. We will then begin with the case in which we assume that the NN direction are 
open and the DD ones are compactified before. The density of states in this case is 

1 I 2 I n=oddM kw \ 10 / \ 10 / 

(57) 

where the string masses now are related to the winding number only 

DD it 2 

M 2 kw = J2^R 2 J+ 2- k. (58) 

3=1 01 

With this density of state is very easy to study the case in which the single string cannot 
probe some DD direction. It is simply accounted by taking some of the winding radii 
to infinity, another important feature is that the IR validity of it allows us to study the 
multiple string gas and then we would be able to analyze the possibility of defining some 
effective dimensions and the existence of energy density scales that probe extra large 
compactified directions. 

Finally we are able to obtain the density of states, for the a null number of non-compact 
NN dimensions. It corresponds to 

= t4 E t [e - (59) 

^2' n=oddM knw \ ) 

The first noticeable feature of this result is that we do not have any continuous depen- 
dence on the energy. The density of states above strictly corresponds to the degeneracy 
of each mass level without any kinematical contribution. On the other hand we know 
that by T-dualizing the NN directions we can describe the system of wrapped D-branes 
as a system of D-particles on a torus. It is interesting to note that if we take the radii of 
the DD direction to infinity we arrive at a system of DO-branes in an already open target 
space. In this case the density of states is reduced to a sum of the mass degeneracy times 
the Dirac-delta function located at each mass level. 

Two relevant limits of this result should be analyzed. If we take the energy to be large 
one we should arrive at a result related to the exponential growth of the degeneracy of 
string mass levels. Doing it in ([59]) we in fact arrive at 

fiW (E = m)~ e p « m m- 9 ' 2 n (m - JJ^J , (60) 

that corresponds to the leading term of (RBI) and to the critical behavior studied in the 
early works on the subject |], §■ As we mention before, because of the absence of 
any kinematical degrees of freedom, the energy here is exactly the mass of the strings. 
What is remarkable is that once again the cut-off establish the range of validity of the 
approximation. 
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On the other hand it is interesting to relate the previous result to the SYM description 
of the D-brane thermodynamics |I| [11], [TJ| and more concretely to the finite temperature 
behavior of Matrix model of M-theory [|T9| , |50|. Taking the massless limit of the open 



string theory on the brane we arrive to the SYM description of the system. On the 
microcanonical side this is implemented by simply taking into account the first term of 
the density of states above. It finally gives a constant density of states that equals the 
polarization of the massless string modes times a delta at zero energy. We can obtain the 
same result by firstly starting by the canonical free energy of the SYM +i theory. It reads 

oo 

F((3)~(3- 1 £ rT 1 . (61) 

n=odd 

If we make the inverse Laplace transformation of Z 1 ^) = —j3F(/3), we easily see that the 
single-object density of states Qgy M (E) exactly matches the massless limit of our result 
in &. 

One final comment. The density of states in fl59|) has the divergent T(|) in the pref- 
actor. The value of this function makes the density of states vanish and this needs some 
interpretation. What is happening is that, because of the compactness of the background 
universe, the spectrum of the string mass is discrete. This means that in a given energy 
range there could be isolated accessible states only, in such a way that we cannot define 
a density of states. In our computation we have supposed the spectrum to be continuous 
and we recover a vanishing for the single string. In principle thermodynamics is not 
affected because, in the microcanonical ensemble, the physical magnitudes come from the 
derivatives of the entropy, which is proportional to the logarithm of the density of states, 
so we can anyway neglect any multiplicative constant. 



6 The number of strings and the entropy. 

In |8|, ||, H 0, [L4| a picture was studied for the phase transition occurring at the Hage- 
dorn temperature. The transition signal was the appearance of a fat string that breaks 
equipartition and absorbs enough energy to excite high mass levels while all the others 
remain massless and form a kind of background, that we called 'the sea' || |§. We would 
like now to see how the number of strings evolves and how the sea and the fat string share 
the energy as the volume varies. To study that, we approximate the high-energy density 
of states by 



with 



n mgh (E, E sea , V) = V(E - E sea )-^e^ E - E ^n sea {E sea ) (62) 

sea - T(dN) N\ { ;{V) ) sea [b6) 

where a is a constant integer that depends on the model and the D-brane background, 
d is the number of open dimensions (testable and continuous at low energies), N is the 
number of strings and E sea is the mean energy of the sea. We can use concepts like average 
energy in spite of working with the microcanonical ensemble because we have separated 
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our microcanonical system into two parts. Each possesses a thermodynamical (physically 
infinite) number of degrees of freedom so that they act as an 'Energy reservoir' for the 
other. That is, the sea of massless strings can be thought of as a gas in the canonical 
ensemble. Throughout the following calculation we shall ignore the thermal fluctuations, 
so that we shall suppose that the energy or the volume are large enough. 

Firstly, we calculate how the energy is shared between both subsystems. The most 
probable distribution is given by the following equation: 



dQ 



High 



The solution to it is 



E, 







E,V 



dN 



(64) 



(65) 



each light string having an average energy d/fin, which is the energy at which the massless 
fields reach the Hagedorn temperature in d spatial dimensions. This result was already 
obtained in || using a model of kinetic theory. 

This is part of the equation of state and it allows us to write 



^■High(E, E sea , V) = Q Hi g h (E, N, V). 



(66) 



This form of the density of states permits to calculate the most probable average 
number of strings for given values of the total energy and the volume and therefore, how 
the energy is shared between the two subsystems. The equation is: 



E 



fat 



dVtuigh 
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dN 


E,V 
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2f3 H . 



-1 - log 



T(d)((d + l)(d- 
2 d - 1 n d / 2 N 



-f(v) 



(67) 



For certain values of V and iV the energy of the fat string becomes small or even 
negative, that is, it ceases to exist. This gives a critical value for /(V) that separates two 
phases with and without the fat string. The phase without fat string is easy to analyze 
because the main contribution comes from the massless string gas. Its density of states is 

N 



n(E,v)=J2 



C(d + i)r(d)(d-r 



N rpdN — 1 



^ T(dN)N\ V 
For large volumes and energies, it yields 



2<i-l 7r £2/2 



V iy E 



(68) 
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d 



d 



2 d-l 7T d/2 



C(d + i)r(d)(d-r 



(69) 



where 



E 



P 



<N> 
<N> 



V 



(70) 
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The subsequent analysis of this expressions depends on the particular form of /(V). 
In the case of closed strings it is V + 1/V, that tends to V away from the selfdual volume. 
Besides, d equals 9 and a equals 1/2. As a result, equation fl6~T|) gives a critical value 
of Peru — 159 strings per a /9//2 . If we now take the thermodynamical limit (V and E 
— > oo) N dynamically diverges maintaining the densities finite. As both densities, p and 
e, grow together, any of them can be taken as the only parameter defining the system. 
The physical picture is as follows: If the density is lower than the critical one, the system 
behaves like a massless gas in ten dimensions and the temperature grows linearly with 
e. When the density reaches the critical value a phase transition occurs and, as the 
correlation length tends to infinity, the system tends to concentrate the energy on one 
fat string with infinite mass. Beyond the transition point, the energy of the fat string 
becomes finite and decreasing with the density. The finiteness of Ef at drive us to formally 
say that its 'Energy density' is zero. We use the inverted commas because one can hardly 
talk about a density for a single object but this is quite a peculiar one because it has 
influence over global properties. Its role is to regulate the temperature so that the bath 
always remains at T#. The system prefers to create more and more light strings rather 
than exciting higher massive modes of the already existing ones. 

The evolution of the temperature is singular. At low energies, the temperature grows 
with increasing specific heat. At the critical point, the temperature is Th and remains so 
thereafter with infinite specific heat. This is a confirmation of the image depicted in ||. 
Fig. (5) shows the qualitative behavior. 

The case of open strings attached to D-branes is quite different because 

f(V) = p-. (71) 

VDD 

The transition depends on how we take to infinite both volumes. If 

> oo (72) 

VDD 

the analysis is identical to that of closed strings, but defining a new density 

P = -TP (73) 

VNN 

It is the density of strings one would define in the T-dual theory where all the directions 
have Neumann boundary conditions. Other possible limit is the opposite one 

VnN - (74) 



V, 



DD 



which is equivalent to taking p to infinite. The energy of the fat string tends to zero but 
is still acts as a regulator so that there is an equilibrium between an unstable fat string 
and the sea that fills the space at Th- 

There is one exceptional case when a = 0, that is, for open strings in a completely 
testable universe. In that case, the breaking of equipartition is quite peculiar. Compare 
figs. (3) and (4). There is never a sharp peak in the density of states of several strings 
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Probability 




Energy 

Figure 3: Probability of the energy distribution for 2-string system in the v = 6 case, at high 
energies. 



signaling the most probable appearance of a fat one that absorbs much more energy than 
the others. Rather, there is a spread in the probability that means that the fluctuations 
in the energy of a single component are comparable to the total energy. This invalidates 
the physical image of the single string dominance of thermodynamics and even prevents 
the correction of it due to what we have called the sea. The correct physical image is a 
gas where equipartition more or less holds but with very large fluctuations in the energy 
of particular strings. This makes the need for the exact calculation stronger because the 
probability of finding a string in virtually any point of the phase space is large. The only 
prediction that can be obtained with the asymptotic expressions is that the asymptotic 
temperature is Hagedorn and that the specific heat is always positive. 

This is the only case in which there is not a phase transition if we take the appearance 
of the fat string as the order parameter. However, at energy densities beyond the Planck 
scale, the system acquires a rather anomalous behavior. For the moment, we cannot 
precise if it corresponds to a higher order phase transition or simply a soft change of 
behavior. 

The difference between the usual microcanonical point of view and this one is the 
fact that here we are letting the number of strings dynamically vary, that is, we are 
working in the (E, V, /i)-ensemble instead of the pure microcanonical. The former is 
more adapted to the string gases because they naturally have a fixed and null chemical 
potential. Somehow, it is a way to gather the different descriptions in terms of finite- 
N gases together. As explained in sections 4 and 5, one can make convolutions of the 
single string density of states and obtain the behavior of different subsystems with a fixed 
number of constituents. The (E, V, /z)-ensemble shows us how the transition between 
them happen. 

The inclusion of the sea as an important contribution to thermodynamics has conse- 
quences on the estimation of the entropy. One has to add the logarithm of (|62|) , that is, 
in the large- N limit 

fVE 9 

S(V, E, N) = logQ(V, E, N) = ION + Nlog-^j^ (75) 



where C is a constant. Substituting the value of N that solves the equation of state (69) 
it is 

10 „ 10 /e^ 9 



S(V,E) = 10N= Y pV = -^-) V (76) 



This is the leading contribution except when a = 0. 
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Let us now write about the pressure. The density of states that we have obtained in 
the previous section has the functional form: 



Q(V,E) = f(E)xg(V) 



(77) 



where 



g(v) 



NN 



V, 



DD 



(78) 



This means that the dependence of the pressure on the temperature is the same as 
that of the free energy. However, its dependence on the volume is not trivial at all. Let 
us recall that, in traditional theories 



P = - 



dF 
dV 



(79) 



However, we know that T-duality is an exact symmetry in type II String Theory. This 
means §J that every physical (measurable) magnitude must be invariant under T-duality. 
This, of course, includes the pressure. The usual definition, that we have written just 
above, is not T-selfdual and therefore it must be corrected and substituted with another. 
We shall use the proposal in |§. 
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This case is even more complicated because the function we are interested in is not 
homogeneous in the volume, so the derivative must be decomposed as 
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with this, and if we assume that all Ri are larger than a', 
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It is interesting to see what happens when we make T-duality act on the free 
so that, for example, all the directions are Neumann, the result would be 
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Figure 4: Probability of the energy distribution for 2-string system in the v = case, at high 
energies. 

where NN and DD are the original Neumann and Dirichlet dimensions. The result is 
similar up to a multiplicative constant term. In fact one can see that what really coincides 
in both points of view of the same theory is the product 

PV oc (NN - DD)^>L (85) 
Vdd 

One remarkable characteristic of the pressure we have calculated is that, in certain 
cases, it vanishes. This happens whenever the number of directions with Neumann and 
Dirichlet boundary conditions coincide. This is only possible if an odd number of dimen- 
sions is untestable. 



7 Conclusions 

In this paper we have studied the thermal behavior of string systems in D-brane back- 
grounds in the microcanonical ensemble. The goal of this kind of analysis is to delucidate 
what is the role of the Hagedorn temperature in the thermodynamics of Dirichlet open 
string theories. As in the closed string cases the necessity of the microcanonical descrip- 
tion comes from the breaking of the canonical ensemble at Hagedorn temperatures. 

Our strategy has been as follows, and is along the lines of ||, |9| . We have begun by a 
computation of the asymptotic density of states of the single string. These computation 
are analogous to those of H, RJ for the closed-string gas and have been recently used for 



the D-brane cases in |18| . We have shown that this procedure corresponds to including an 
IR cut-off A 2 that eliminates the contribution of the massless string modes. Consequently 
the densities of states obtained are not valid for the low energetic regime. We have studied 
the range of validity of this approximation in terms of the energy of the system, the radii 
of the compact dimensions and the value of A 2 . The result is a kind of phase space where 
different parts are best described by particular approximations. Only two of these parts 
are well described by the asymptotic density of states, approximated by the Hagedorn 
term only. Passing from one of them to the other the number of effective dimensions of 
the target space the string gas lives on changes. We claimed that in order to completely 
know the thermal behavior it is necessary to use a complete, IR regular density of states. 

On the other hand the thermodynamics of the string gas is not exactly described 
by the single string system, so we need to study the multiple string density of states. 
We should follow the convolution prescription of pi. However it is not possible to do it 
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using the high energy densities of states of the single string because of their irregular IR 
behavior. Finally one is lead to look for a complete density of states. 

Before computing the desired density of states we have analyzed the thermal properties 
of the single string. We have shown that there is a negative specific heat phase in all cases 
except for the D + v = 9 one. In the latter we have an IR definite phase but we think that 
at these energies the density of states does not include all the information. It is in fact 
obvious that the corrections to the leading behavior are strongly dependent of A 2 , and 
that way nothing warrants that we are including all the IR degrees of freedom of the 
theory. The only extractable conclusion is that after including all the IR information we 
will still have a positive specific heat phase, and Hagedorn will be a asymptotic maximum 
temperature of the system. 

We have calculated an exact expression for the single string density of states in section 
5. It is got in the form of a series with a finite number of terms for finite values of the 
energy. The results are analogous of those in 0,0. With those instruments, we are 
able to study the thermodynamics of the systems and also to match our results with 
those obtained using the SYM p+ i description of the D-branes dynamics. We have used 
the D + v = 9 cases as some special examples. In one of them, DD = 9, we are able 
to connect with the Matrix model microcanonical description, in the Born-Oppenheimer 
approximation. 

In particular, we study the distribution of energy among the constituents. This gives 
us physical images that help us test models that had been proposed in earlier works 
f§, |9], ||. When a fat string with more energy than the rest appears and, therefore, 
equipartition is broken, a phase transition occurs. We have been able to distinguish the 
cases in which this happens and in which it does not. We have seen how this affects the 
temperature and the specific heat value, the case in which the phase transition occurs has 
been plotted in Fig. (5). When several of this systems are put in thermal contact, that 
is, when open strings attached to different D-branes and free closed strings thermally 
interact, the behavior of the global system is dominated by the most degenerate phase. 
This is the one that gives lower temperature and smaller specific heat for a value of the 



energy: open strings with D + v = 9. A similar picture has already been obtained in [18 



the difference between our picture and theirs is that in ours all subsystems have T# as the 
maximum temperature. There is not any 'non-limiting' system with temperatures higher 
that Th- Our result is that systems can be classified between those with phase transition 
and those without it. The systems without phase transition are the closed strings in a 
completely compactified space || and open strings with D + v = 9. Whenever there is 
at least one open direction, the only system that does not suffer phase transition is that 
of open strings with D + v = 9. When one considers the whole system made up with 
closed strings and open strings of various types (various values of u,DD and NN), the 
most degenerate system dominates the thermodynamics. It is always the one without 
transition. This means that the existence of one D-brane prevents any phase transition, 
except when some DD are open. 

Along all these considerations we have considered all the untestable DD direction to 
be already open. If we had not done it we would have found some of the previous systems 
in thermal contact having the phase transition. It corresponds to systems that stand in 
region III of fig. (2). This systems have the energy per string (p N ) constraint by 
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Figure 5: Temperature versus energy for the multiple string systems that reach Hagedorn 
temperature through the long string nucleation 

On the other hand we know that the phase transition occurs at E/N ~ DTh = pu- 
Two possibilities come. If pu < Pu the systems suffers a phase transition. Else we need 
to include all winding states and finally we pass to the region / where we do not find 
any phase trasition. We think that in the thermodynamic limit, where p^ goes to zero, a 
phase transition is prevented. 

Regarding the pressure, we calculate it in several cases and see the peculiarities of 
string systems using the String-Theory definition of it proposed in ||. 

Finally, we use the 'fat string-sea' model to make a calculation in a genuinely ther- 
modynamic limit, and see the evolution of the temperature as a function of string and 
energy densities instead of absolute magnitudes. 
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